Melting of 'porous' vortex matter 
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Bitter decoration and magneto-optical studies reveal that in heavy-ion irradiated superconduc- 
tors, a 'porous' vortex matter is formed when vortices outnumber columnar defects (CDs). In this 
state ordered vortex crystallites are embedded in the 'pores' of a rigid matrix of vortices pinned 
on CDs. The crystallites melt through a first-order transition while the matrix remains solid. The 
melting temperature increases with density of CDs and eventually turns into a continuous transition. 
At high temperatures a sharp kink in the melting line is found, signaling an abrupt change from 
crystallite melting to melting of the rigid matrix. 

PACS numbers: 74.60.Ec, 74.60.Ge, 74.80.-g, 74.72.Hs 



Melting of heterogeneous systems, and in particular 
of nanocrystals embedded in porous rigid matrices, is 
a complex process with many uncontrolled parameters. 
Metal and semiconductor nanocrystals with free surfaces, 
for example, usually show a decrease in their melting 
temperature with decreasing size [1], whereas nanocrys- 
tals encapsulated in a porous matrix often display an 
increase in melting temperature [2]. Although the con- 
tribution of the different factors is still a matter of de- 
bate, the melting process is known to depend on the size, 
dimensionality, material properties of the nanocrystals 
and the matrix, as well as the interface energies between 
the materials [1,2]. In this work we investigate an analo- 
gous, but a more controllable composite system, which is 
a 'porous' vortex matter consisting of vortex nanocrys- 
tals encapsulated in a matrix of strongly pinned vortices. 
As shown below, this system is present in the commonly 
heavy-ion irradiated superconductors when the vortices 
outnumber the columnar defects (CDs). The rigid ma- 
trix is created by vortices localized on the network of 
random CDs, while the softer nanocrystals are formed 
within the 'pores' of this matrix by the interstitial vor- 
tices. The size of the nanocrystals can be readily varied 
from several hundred down to a few vortices by changing 
the applied field or the density of CDs. We find that this 
composite vortex matter reveals a number of intriguing 
mechanisms: Similarly to the metallic nanocrystals in a 
matrix, we observe for the first time a pronounced up- 
ward shift in the vortex melting temperature T m , while 
preserving the first-order nature of the transition (FOT) . 
With increasing density of CDs, the size of the pores de- 
creases, resulting in a larger shift in T m . We also find a 
critical point at which the FOT changes into a contin- 
uous melting. Moreover, the crystallites can melt while 
the matrix remains rigid. As a result, at high tempera- 
tures we find an abrupt breakdown in the upward shift 



of T m and a sharp kink in the FOT line, which appar- 
ently result from the collapse of the matrix due to vortex 
depinning from the CDs. 

The reported findings were obtained using Bitter dec- 
oration and differential magneto-optical (MO) [3] tech- 
niques. High quality Bi 2 Sr 2 CaCu208 (BSCCO) crystals 
(T c « 89 K) were covered by various patterned masks 
and irradiated at GANIL by 1 GeV Pb ions with doses 
corresponding to matching fields of B$ — 5, 10, 20, and 
50 G. Figure la shows schematically one of these masks 
which results in the formation of CDs in BSCCO crystals 
only within the circular apertures of about 90 fim diam- 
eter. This patterning allows very sensitive simultaneous 
comparison of the vortex structure and the local melt- 
ing processes in adjacent irradiated and pristine regions, 
which is not possible by other methods. 

In the absence of CDs the vortices form the Bragg glass 
phase [4] which has quasi-long-range order with no topo- 
logical defects, as seen in the pristine part of the mag- 
netic decoration image in Figs, lb and lc. What hap- 
pens to this phase when sparse CDs are added? The 
irradiated region in Fig. lb shows that it is no longer 
Bragg glass since it has significant amount of topologi- 
cal defects (solid circles in Fig. lc) and no orientational 
long-range order [5]. It is also not an amorphous or 
glass phase in the usual sense, nor is it a simple poly- 
crystal as discussed below. In the presence of CDs the 
Bose glass (BG) theory [6-8] is usually applied, which 
describes the vortex matter in terms of anisotropic ho- 
mogeneously pinned medium. Such a description is ade- 
quate for the common situation of high irradiation doses 
(Bcf, > B), in which all the vortices reside on CDs and 
the vortex pinning energies are comparable. In contrast, 
we investigate here mainly the opposite extreme of low 
doses in which the vortices greatly outnumber the CDs, 
i.e., B 3> Bcj,. In this case the system is inherently hetero- 
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geneous, consisting of two vortex populations with well 
separated characteristic energies: The vortices residing 
on CDs are strongly pinned and form a rigid network 
or matrix, whereas the interstitial vortices are localized 
by significantly weaker elastic interactions and form rel- 
atively soft crystallites within the pores of the matrix. 
In the following we refer to this state (upper parts of 
Figs, lb and lc) as porous vortex matter in order to em- 
phasize the important consequences of the heterogeneous 
structure. 

Figure 2 shows several frames from a 'movie' [9] of the 
melting process as a function of temperature T at two 
fields, 30 and 60 G. Each frame is obtained by taking 
the difference between the MO images at T + 0.15 K and 
T — 0.15 K and averaging a large number of such dif- 
ferential images, as described previously [3]. The bright 
features show the regions in the sample that undergo a 
FOT within the temperature interval of 0.3 K at the in- 
dicated T, and the intensity of this bright paramagnetic 
signal is the equilibrium magnetization step AB at the 
transition [3,10]. Figure 2a shows the nucleation of the 
liquid phase in the form of bright inclined strips in the 
central pristine region which arise from intrinsic sample 
disorder (see Rcf. [3]). With increasing T (Fig. 2b) the 
liquid expands, remarkably avoiding the irradiated aper- 
tures. In Fig. 2c the entire central pristine part of the 
sample is liquid, while the apertures with = 20 G are 
still solid. In Fig. 2d the central apertures melt at 82.45 
K, which is about 1 K above T m of the adjacent surround- 
ing pristine regions in Fig. 2b. The apertures closer to 
the sample edges begin to melt in Fig. 2e. The FOT is 
equally strong in the irradiated and pristine regions: The 
AB step derived from the paramagnetic melting signal 
[10] is the same in Figs. 2d and 2b. Also the width of the 
local melting transition is the same within our resolution, 
i.e., each point in the sample melts within 0.3 K or less. 
This is the first direct observation of an upward shift of 
the FOT by correlated disorder. Note that for small B^ 
the melting in the irradiated apertures occurs while the 
melting in remote pristine regions is still in progress, and 
therefore the quantitative measurement of the shift in T m 
has to be determined by comparing neighboring irradi- 
ated and pristine regions and cannot be readily detected 
by global techniques. 

The melting process at 60 G (Fig. 2, second row) re- 
veals two important differences. First, the shift of the 
melting temperature, AT mi is about 4 K (difference be- 
tween Figs. 2i and 2f) which is much larger as compared 
to 1 K at 30 G. In Fig. 2h, for example, the entire pristine 
sample has melted while the irradiated apertures are still 
solid. Second, the brightness of the paramagnetic melt- 
ing signal AB in the apertures in Fig. 2i is much lower 
than in the pristine sample, and moreover, the melting 
in each aperture is broadened over several frames, as ex- 
emplified in Fig. 2j. At still higher fields, above 100 G, 
no paramagnetic FOT signal is detected in the irradiated 



apertures with B^ — 20 G. 

In order to investigate the phase transition over a wider 
range of fields we have used differential MO imaging in 
which the field is modulated by 1 G (Fig. 3) instead 
of the T modulation. In addition to detecting the FOT 
[3] , this method provides a very sensitive measurement of 
the irreversibility line at a very low effective frequency of 
about 0.1 Hz. Figure 3 demonstrates the determination 
of the irreversibility line in the region where no FOT is 
present. In Fig. 3a the entire pristine sample is in the 
liquid state, while the apertures are still solid. The aper- 
tures appear black, showing that the external field mod- 
ulation is shielded due to the enhanced pinning. Upon 
increasing the field the black apertures disappear (Figs. 
3b and 3c) revealing the value of the local irreversibility 
field. 

Figure 4 shows the location of the onset of the FOT 
for B<j, = 5, 10, 20, and 50 G obtained by T modulation 
(solid symbols) and of the irreversibility line obtained 
by field modulation (open symbols). The solid lines are 
guides to the eye for the FOT lines which terminate at 
the corresponding critical points. The irreversibility data 
coincide with the FOT line below the critical point and 
smoothly extrapolate the location of the transition line 
to higher fields. The first interesting observation here is 
that although the structures of the porous vortex matter 
and of the Bragg glass are very different (Fig. 1), their 
phase diagrams for B^ = 5 G in Fig. 4 are almost identi- 
cal. The melting remains FOT in most of the field range 
and the melting line is shifted only slightly. This brings 
us to an important conclusion that the quasi-long-range 
order that characterizes the Bragg glass is not an essen- 
tial requirement for the existence of a FOT [11], and the 
presence of the short-range order within the crystallites 
is apparently sufficient. 

We can understand the upward shift in B m {T) in Fig. 
4 by generalizing the concept of the cage model. In a 
pure system each vortex is confined in a potential cage 
arising from the elastic interactions with its neighbors. 
Melting occurs when the transverse thermal fluctuations 
of vortices (ut) reach a certain fraction cl of the lattice 
spacing ao. The pores of the matrix are vertical cylinders 
which provide an additional confining cage potential to 
the interstitial vortices. This enhanced rigidity reduces 
(ut), and hence stabilizes the solid crystallites within 
the pores, thereby shifting B m (T) upwards. The shift in 
B m (T) grows with B^ due to the decrease in the size of 
the pores. Since the pinning energy of the CDs is usually 
significantly larger than the elastic energy, the melting of 
the crystallites may occur without the destruction of the 
matrix, which remains solid up to a higher temperature 
as described below. The rigid matrix may thus coexist 
with an interstitial liquid as shown theoretically [8] and 
in numerical simulations [12]. 

Upon increasing the density of CDs the FOT is weak- 
ened and it eventually transforms into a continuous tran- 
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sition in Fig. 4. This transformation is seen more clearly 
by varying the temperature along a given B m (T) line. 
The inset to Fig. 4, left axis, shows the height of the 
equilibrium magnetization step AB vs. T for = 20 
G. As the temperature is decreased AB drops and even- 
tually vanishes at a critical point (CP) below which no 
discontinuity is found and the melting becomes contin- 
uous. In addition, the width of the FOT, ST m , shows 
a significant broadening on approaching the CP (Fig. 4 
inset, right axis). This indicates that it is a true CP at 
which the FOT transforms into a continuous transition 
as expected theoretically [13]. This CP is unlike the ap- 
parent point-disorder-induced CP in BSCCO at lower T 
where no appreciable broadening was observed [14] , and 
where the FOT was found to be obscured by pinning- 
induced irreversibility in the liquid phase [15]. Here, in 
contrast, the liquid is fully reversible and no such exper- 
imental limitations exist. Figure 4 shows that the CP is 
shifted to lower fields with B^ [16]. It is interesting to 
note that in YBCO the CP was observed to shift upward 
with CDs [17]. There, however, the CP is induced by 
the point disorder in the pristine crystals and the CDs 
just modify slightly its position by reducing the point- 
disorder-induced vortex meandering. In our case, in con- 
trast, the CDs apparently create a new critical point that 
does not exist in pristine crystals. 

We now discuss another key finding which is a kink in 
the B m (T) lines marked by the arrows in Fig. 4. This 
kink becomes very prominent when the temperature shift 
AT m between the irradiated and pristine melting lines 
is plotted as shown in Fig. 5. The inset displays the 
corresponding upward shift in field AB m . In BG the- 
ory a kink in the transition line is expected to occur at 
Bk ~ B^ [7,8]. We argue that the kink in Fig. 5 is a new 
feature of completely different nature, which reflects the 
collapse of the rigid matrix. The BG kink is experimen- 
tally found to occur at Bk/B^ of 1/6 to 1 for large B^, 
and Bk usually scales linearly with B^ [18,19]. The kink 
in Fig. 5 occurs in the opposite regime of Bk/B^ m 2 to 
8, and it does not scale with B^. At high B^ the BG kink 
is a broad feature that was found to occur only along a 
glass transition line [18,19]. In contrast, here the kink is 
extremely sharp, and it is the first observation of a kink 
that occurs along a FOT line. Finally, the BG kink is a 
transition from a strong influence of the CDs below the 
kink to a much weaker effect of CDs above the kink [7,8] . 
Here, in contrast, the situation is just the opposite, as 
seen in the inset to Fig. 5. The effect of CDs is large 
at fields above the kink, where the shift AB m is almost 
constant and approximately equal to B^, and it collapses 
rapidly at T > Tk- This collapse occurs at Tk — 75 K, 
with little dependence on B^. It is also interesting to note 
that for B^ < 20 G the kink occurs along the FOT line, 
while for B^ = 50 G it falls in the region of continuous 
transition. 

In Fig. 4 three regions are enumerated with respect 



to B^ = 50 G curve. In region 1 the crystallites in the 
pores are stabilized by the rigid matrix and melt at B!™ r , 
well above the pristine B m (T). In region 2 the pores are 
liquid while the matrix remains intact, and hence B]™ r 
reflects a melting transition of the softer of the two sub- 
stances in a heterogeneous medium. This unconventional 
melting cannot be described by a single set of parameters 
since it depends on the properties of both the pores and 
the matrix, each described by a different set of param- 
eters with different field and temperature dependencies. 
As a result, B^ r (T) does not extrapolate to T c as com- 
monly expected, but rather well above it. To the best of 
our knowledge this is the first observation of a melting 
or irreversibility line that has such an unusual property 
at high temperatures. This heterogeneous process col- 
lapses sharply at Tk, above which the matrix rapidly de- 
localizes, resulting in a homogeneous liquid in region 3. 
The B™ fa line thus describes the melting process of the 
matrix due to dcpinning from the CDs, which leads to 
immediate melting of the superheated crystallites. Thus 
in contrast to the common homogeneous BG description, 
the Bf^ r (T) and B™ tx (T) lines originate from two inde- 
pendent heterogeneous processes, melting of crystallites 
within a rigid matrix and the collapse of the matrix it- 
self, resulting in a sharp kink at the intersection point. 
We expect the B™ tx (T) line to extend also into the liq- 
uid phase [8], like the dashed line in Fig. 4, separating 
the interstitial liquid within a solid matrix in region 2 
from the homogeneous liquid 3. Our MO measurements 
cannot detect this line since both these regions are fully 
reversible. In region 2, however, we expect a higher c-axis 
correlation than in region 3. Recent Josephson plasma 
resonance studies of BSCCO with low B$ indeed find 
a recouping transition within the liquid phase at which 
an enhancement in the interlaycr phase coherence is ob- 
served [19]. This extension of the B™ tx line should also 
be detectable by transport measurements since the flux- 
flow resistance of the interstitial liquid should be lower 
than of the homogeneous liquid [8] . These studies will be 
the subject of future work. 

In summary, when vortices outnumber CDs, hetero- 
geneity rather than the average properties of the lattice 
has to be taken into account for a proper description 
of the structure and the thermodynamic behavior of the 
vortex matter. We have found evidence for two mech- 
anisms: melting of superheated crystallites within the 
pores of a solid matrix and the destruction of the rigid 
matrix. The intersection point of these two independent 
processes results in a sharp kink in the observed melting 
line. The heterogeneous melting can be either first-order 
or continuous, depending on temperature and the density 
of CDs. The porous vortex matter may thus provide a 
tunable model system for general comprehension of melt- 
ing of nanocrystals in porous materials. 
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FIGURE CAPTIONS 

Fig. 1. (a) Schematic of BSCCO crystal irradiated through a mask with an array of circular apertures, (b) Bitter 
decoration image (B = 40 G, T = 4.2 K) showing pristine region at the bottom (Bragg glass) and a section of the 
irradiated aperture (B^ = 10 G) on top where porous vortex matter consisting of ordered crystallites embedded in a 
rigid matrix is formed, (c) Corresponding locations of six-fold (open circles) and of five- and seven-fold coordinated 
vortices (solid circles) obtained by Delaunay triangulation. 

Fig. 2. Melting process as a function of T at fields of 30 and 60 G in BSCCO crystal 0.65 x 0.45 x 0.01 mm 3 using 
differential MO imaging with T modulation of 0.3 K. The bright regions are the areas that undergo melting within 
the 0.3 K interval, (a) — (c) and (f) — (h): melting of the pristine regions while the irradiated apertures are still 
solid, (d), (e), (i), and (j): melting of the irradiated apertures {B$ = 20 G). The area outside the crystal is blackened. 
Color movies of the melting process are available at http:/ /www. weizmann.ac.il/home/fnsup/. 

Fig. 3. Differential MO images using field modulation of 1 G in BSCCO crystal of Fig. 2 at three fields at T = 68.4 
K. (a) All the pristine regions are liquid while the irradiated apertures {B^ = 20 G) are still solid and irreversible, 
and hence appear black due to shielding, (b) Partial reversibility of the central apertures, (c) Central apertures are 
liquid while the apertures closer to the edges begin to melt. 

Fig. 4. The melting lines B m (T) of the pristine and irradiated regions with indicated B^. Solid (open) symbols are 
temperature (field) modulation data showing the location of the FOT (irreversibility line). Solid (dotted) lines are 
guides to the eye of the first-order (continuous) transitions that terminate at the critical points Q. Inset: The height 
of the FOT equilibrium magnetization step AB which vanishes at the CP, and the local width ST m of the FOT vs. 
T. 

Fig. 5. The shift AT m in the melting temperature for different B^ with respect to the pristine T m at various fields 
B. Inset: The upward shift AB m in the melting field vs. T m . 
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